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ABSTRACT 

This paper is the third one in a series, intended to investigate how the observed kine- 
matics of elliptical galaxies are affected by dust attenuation. In Paper I and Paper II, 
we investigated the effects of dust absorption; here we extend our modelling in order 
to include the effects of scattering. We describe how kinematical information can be 
combined with the radiative transfer equation, and present a Monte Carlo code that 
can handle kinematical information in an elegant way. 

Compared to the case where only absorption is taken into account, we find that 
dust attenuation considerably affects the observed kinematics when scattering is in- 
cluded. For the central lines of sight, dust can either decrease or increase the central 
observed velocity dispersion. The most important effect of dust attenuation, however, 
is found at large projected radii. The kinematics at these lines of sight are strongly 
affected by photons scattered into these lines of sight, which were emitted by high- 
velocity stars in the central regions of the galaxy. These photons bias the LOSVDs 
towards high line-of-sight velocities, and significantly increase the observed velocity 
dispersion and LOSVD shape parameters. These effects are similar to the expected 
kinematical signature of a dark matter halo, such that dust attenuation may form 
an alternative explanation for the usual stellar kinematical evidence for dark matter 
halos around elliptical galaxies. 

We apply our results to discuss several other topics in galactic dynamics, where 
we feel dust attenuation should be taken into account. In particular, we argue that 
the kinematics observed at various wavelengths can help to constrain the spatial dis- 
tribution of dust in elliptical galaxies. 

Key words: dust, extinction - galaxies: elliptical and lenticular, cD - galaxies: 
kinematics and dynamics - radiative transfer - scattering 



1 INTRODUCTION 

It is now generally accepted that elliptical galaxies are com- 
plicated objects, containing a variety of stellar populations, 
with a still poorly constrained dark matter content and a 
multi-component interstellar medium. The latter, and in 
particular the interstellar dust component, is not well un- 
derstood. Nevertheless, knowledge about the dust compo- 
nent in galaxies is of fundamental importance. On the one 
hand, interstellar dust plays an active role in e.g. interstellar 
chemistry and star formation, and is therefore a key ingre- 
dient to understand galaxy structure and evolution. On the 
other hand, dust is very effective in absorbing and scatter- 
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ing UV and optical light, and knowledge about its presence, 
quantity and properties is necessary to correctly interpret 
any observable. 

The far-infrared fluxes detected in the eighties by the 
IRAS satellite at 60 fim and 100 /im unveiled the presence of 
a substantial amount of dust in elliptical galaxies (Jura 1986; 
Bally & Thronson 1989; Knapp, Gunn & Wynn- Williams 
1989). At that time, dust was already observed in the form 
of dust lanes and patches for a number of early-type galaxies 
(Hawarden et al. 1981; Ebneter & Balick 1985; Veron-Cetty 
& Veron 1988). Recent surveys indicate that dust extinction 
features are present in a large fraction of early-type galaxies 
(van Dokkum & Franx 1995; Ferrari et al. 1999; Tomita et 
al. 2000; Tran et al. 2001). These dust features, however, 
can not account for the high IRAS fluxes: the dust masses 
estimated from the IRAS measurements exceed those esti- 
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mated from integrating the extinction features by nearly an 
order of magnitude (Goudfrooij & de Jong 1995). Moreover, 
the IRAS dust mass estimates are a lower limit for the true 
dust masses, because IRAS is not sensitive to cold dust, 
which emits the bulk of its radiation longwards of 100 /im. 
Calculations with a more realistic dust temperature distri- 
bution (Merluzzi 1998) and submillimeter observations (Fich 
& Hodge 1993; Wiklind & Henkel 1995) may indicate dust 
masses up to a magnitude higher than the IRAS estimates. 
The only plausible way to solve this dust mass discrepancy 
is to assume that the interstellar dust in elliptical galaxies 
exists as a two-component medium: the less massive compo- 
nent is optically visible in the form of dust lanes, whereas 
the more massive one is distributed diffusely over the galaxy 
(Goudfrooij & de Jong 1995). 

Such a diffusely distributed dust distribution can be 
traced in two ways. The most obvious way is to map the 
thermal emission of the dust with spatially resolved far- 
infrared or submillimeter observations. A first effort to do 
this was undertaken by Leeuw, Sansom & Robson (2000), 
who used SCUBA imaging at 850 /im to look for thermal 
emission from dust in NGC 4374. Unfortunately, they could 
not detect a diffuse dust component, and found that the 
bulk of the 850 /jm core emission is more likely the result 
of synchrotron radiation. Deeper imaging covering a larger 
wavelength range is necessary to constrain the spatial dis- 
tribution of dust this way. Hopefully, the new generation 
of far-infrared and submillimeter instrumentation such as 
SIRTF and ALMA will help to clarify this issue. 

A second way in which the diffuse dust component in 
elliptical galaxies can be traced is by colour gradients: be- 
cause the extinction efficiency of dust grains decreases with 
wavelength, a diffuse dust component is expected to generate 
bluer colours at larger projected radii. Goudfrooij & de Jong 
(1995) and Wise & Silva (1996) found that the dust distri- 
butions, necessary to create the colour gradients observed in 
a sample of elliptical galaxies, are in reasonable agreement 
with their observed integrated FIR fluxes. However, various 
arguments seriously complicate the interpretation of colour 
gradients. Foremost, dust attenuation^] is not the only pro- 
cess that generates colour gradients: also age and metallicity 
variations can cause broadband colour gradients. Recently, 
Michard (2000) argued that the mean observed colour gra- 
dient ratios of a sample of elliptical galaxies are more likely 
the result of metallicity than of diffuse dust. And even if 
the dust attenuation were the only process responsible for 
the generation of colour gradients, tracing and quantifying 
the amount of dust would still be complicated. First, the 
dust is not located between the source and the observer, 
but well mixed with the stars. As a result, the amount of 
dust is not simply proportional to the amount of reddening 
(Disney, Davies & Phillipps 1989). Second, bluing due to 
scattering partly compensates the effects of reddening due 
to absorption, which suppresses the formation of large broad 
band colour gradients even if a substantial amount of dust 
is present (Witt, Thronson & Capuano 1992). 

We have set up a program to investigate the effects of 
diffuse dust on the observed kinematics of elliptical galax- 

1 We refer to dust attenuation as the combined effect of absorp- 
tion and scattering. 



ies. In the first two papers of this series (Baes & Dejonghe 
2000, hereafter Paper I; Baes, Dejonghe & De Rijcke 2000, 
hereafter Paper II) we investigated how the light profile and 
the observed kinematics of elliptical galaxies are affected 
by dust absorption. In this series we extend our models by 
including the effect of scattering. In Section ^| we explain 
how the processes of dust absorption and scattering can af- 
fect the observables of galaxies. In particular, we explain 
how kinematical information can be included in the radia- 
tive transfer equation. In Section [], we argue that a Monte 
Carlo method is the most straightforward way to do this. 
We describe the code we developed, with a special emphasis 
on the calculation of the observed kinematics. We use this 
code to investigate the effects of dust attenuation on the 
observed kinematics in elliptical galaxies. In Section ^| we 
present a set of simple elliptical galaxy models, consisting of 
a stellar and a dust component. We demonstrate how dust 
attenuation affects the light profile and observed kinemat- 
ics of these models in Section ^. Section |^ is devoted to a 
discussion of our results, and finally, Section H sums up. 



2 RADIATIVE TRANSFER IN DUSTY 
GALAXIES 

2.1 The general radiative transfer equation 

The basis for any study of attenuation is the radiative trans- 
fer equation (RTE) , which statistically describes the interac- 
tion between matter and light. In a general form, the time- 
independent RTE can be written as (Chandrasekhar 1960; 
Mihalas 1978) 

^(r,fc)=i(r,fc)-K(r)I(r,fc), (1) 

where s is the path length, and I(r,k) represents the inten- 
sity at a position r into a direction fe. The right-hand side 
of the equation contains two terms that represent how the 
radiation field changes as a result of interactions with the 
matter. The first term, the total emissivity j(r, fe) accounts 
for the sources of the radiation field. The second term, where 
Ac(r) is the total opacity coefficient, represents the sinks. The 
specific form of both terms depends on which physical pro- 
cesses are taken into account, (e.g. stellar emission, absorp- 
tion, scattering) and they can be complicated functions that 
depend on the intensity itself. 

In a galaxy without dust attenuation, there are no sink 
terms, and the only source in the radiation field is the 
(isotropic) emission of photons by stars, such that j(r, fe) = 
l(r) where £(r) represents the stellar emissivity. If we take 
absorption by dust grains into account, a sink term must be 
added, which accounts for the loss of photons from the ra- 
diation field.0 The total opacity then equals the absorption 
coefficient, i.e. «(r) = « a b s (r). In either of these situations, 
the RTE is a ordinary differential equation, and it can be 
solved by an integration along the line of sight. 

2 In fact, dust absorption will also account for an additional 
source term, because the energy absorbed by grains will be re- 
emitted at infrared wavelengths. Because we are primarily in- 
terested in the optical and near-infrared regimes, however, this 
source term can safely be neglected. 
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Dust grains do not only absorb photons, they also scat- 
ter them, i.e. a number of photons are, as a result of an 
interaction with a dust grain, removed from their path and 
sent into another direction. This physical process is not a 
rare phenomenon: for typical Milky Way dust grains, the 
probability for scattering even slightly exceeds the probabil- 
ity for absorption in the optical wavelength range (Table g). 
It is therefore obvious that scattering should be included in 
radiative transfer calculations. This will add two extra terms 
to the RTE. The first one is a sink term that accounts for the 
loss of photons scattered out of the beam. The total opacity 
will hence be the sum of the absorption coefficient and the 
scattering coefficient, n(r) = K a b s (r) + n sca _(r). The second 
extra term is a source term that characterizes the gain of 
photons scattered into the beam. More precisely, this term 
will contain the contribution of photons that had another di- 
rection fe', but are now scattered into the direction k. Gener- 
ally, the distribution of angles after a scattering process are 
described by a scattering phase function <fr(r,k,k'), which 
describes the probability that a photon which comes from 
the direction k! and is scattered at r, will have k as its new 
direction. By convention, it is normalized as 



// 



for all r and all k. 



(2) 



The fraction of the intensity scattered from a solid angle 
dfi' around an arbitrary direction k' into the direction k 
will therefore equal K sca .(r) I(r, k') $(r, k, k') dfi'/47r, such 
that the extra source term that has to be added to the RTE 
reads 



I(r,k')<f>(r,k,k') 



(3) 



It is convenient to introduce the scattering albedo ui as the 
ratio of the scattering coefficient to the total opacity coeffi- 
cient, 



ui(r) = r-V^- = 1 - 



«abs(r) 



(4) 



We finally find for the RTE, 
j-(r,k) = l(r) - k(t) I(r,k) 

+ w(r) K (r) JJ^I{r,k')*(r,k,k'). (5) 

The inclusion of scattering turns the RTE into a integro- 
differential equation, far more complicated than the ordi- 
nary differential equation if only absorption is taken into 
account. In particular, the scattering term is responsible for 
the coupling of the RTE along different paths: due to the 
integration over the angle, we cannot solve the RTE for a 
single path, but we have to solve it for all paths at the same 
time. 



2.2 The geometry 

The complexity of the RTE does not only depend on which 
physical processes are taken into account (the right-hand 
side) , but also on the geometry of the system. Indeed, the 
path length appearing in the left-hand side of the RTE is in 
general a function of position and direction, i.e. s = s(r, fe). 
As a consequence, the (time-independent) RTE is a partial 



differential equation with five independent coordinates (see 
e.g. Mihalas 1978). This complexity, however, is reduced if 
the system has symmetries. For example, in an axially sym- 
metric geometry, the azimuthal dependence vanishes, such 
that only four independent coordinates need to be consid- 
ered. 

Particularly interesting is the spherical geometry, be- 
cause in that case only two coordinates remain: the radius 
r and the cosine fi of the angle between the direction fc and 
the local radial direction. The aim of this paper is to pro- 
vide a global picture of the effects of dust attenuation on the 
observed kinematics of elliptical galaxies, rather than mod- 
elling a specific, possibly geometrically complex object. For 
this goal, the assumption of spherical symmetry is satisfac- 
tory. 

There are no ways to solve the general RTE (H) analyti- 
cally, but several techniques have been developed to solve it 
numerically. And because many astrophysical systems can in 
a first approximation be considered spherically symmetric, 
many efforts have yet been spent on the RTE in a spherical 
geometry. The pioneering work started nearly 70 years ago 
(Kosirev 1934; Chandrasekhar 1934), and, in particular since 
the 1970s, a vast number of different approaches has been de- 
veloped (e.g. Hummer & Rybicki 1971; Schmid-Burgk 1975; 
Witt 1977; Flannery, Roberge & Rybicki 1980; Yorke 1980; 
Rowan-Robinson 1980; Rogers & Martin 1984; Peraiah & 
Varghese 1985; Gros, Crivellari & Simonneau 1997). Most 
of these techniques, however, are not suitable for our needs, 
because it seems hard (and at least not obvious) to extend 
them such that they can handle kinematical information. 
Only for Monte Carlo methods, the kinematical information 
could be included in an elegant way. 

2.3 Including kinematical information 

As long as scattering is not taken into account, the inclusion 
of velocity information in the RTE is rather straightforward. 
Indeed, instead of taking into account the entire stellar emis- 
sivity, we can just consider the light from those stars whose 
velocity component in the direction of the observer equals 
v\\ . The emissivity of these stars is given by £(r) 4>(r, fc, v\\), 
where 4>{r, fe, v«) is the spatial LOSVD, i.e. the probability 
for a star at position r to have a velocity component vu in 
the direction fe (Appendix A). Solving the RTE with this 
emissivity, we obtain the LOSVDs observed in the plane 
of the sky. Hence, if either dust attenuation is completely 
neglected or only absorption by dust grains is taken into 
account, the observed LOSVDs can be calculated from the 
spatial LOSVDs through a single integration along the line 
of sight. Similar relations hold between the moments of the 
distribution function, and the moments of the LOSVDs, e.g. 
the observed velocity dispersion profile can be determined 
from the intrinsic velocity dispersions by a simple integra- 
tion along the line of sight. For more details, we refer to 
Paper I. 

When scattering is included, however, the inclusion of 
kinematical information becomes much more complicated, 
because the velocities of both the stars and the dust grains 
that scatter their light need to be taken into account. Con- 
sider, for example, Fig. |l|. A photon with rest wavelength Ao 
is emitted by a star moving with a velocity u* into a direc- 
tion feo (in a reference frame centered on the galaxy center). 
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observer 



Figure 1. The inclusion of kinematical information in the ra- 
diative transfer equation. This figure shows the trajectory of a 
photon through the galaxy (undulating line): after being emitted 
by the star, it is scattered once before it leaves the galaxy. See 
text for more details. 



Assume this photon is scattered by a dust grain moving with 
velocity Vd, and then scattered towards the observer. The 
wavelength "detected" by the dust grain is determined by 
the relative velocity between star and dust grain, i.e. the 
dust grain detects a wavelength 

(w* - v d ) ■ k 



A = A 



(6) 



The dust grains scatters the photon coherently into the di- 
rection fc bs, such that the wavelength detected by the ob- 
server equals 



A = A 



1 - 



(«» - v d ) ■ k 



1 - 



(v d - 



S )-k 



(7) 



where v sys is the system velocity with respect to the ob- 
server. The line-of-sight velocity Vu detected by the observer 
can be found from the relation 

feobs 



A = A 



V\\ + V B 



(8) 



which defines v\\ and already incorporates the system veloc- 
ity. Because the velocities of dust grains and stars in galax- 
ies are very small compared the speed of light, second order 
terms in v/c can be safely neglected, such that with ([?]) and 

(1). 



feo) 



(9) 



W|| = V, ■ feo + V d ■ (fe oba 

When more scattering events are involved in the photon's 
trajectory, the relative velocities of each pair of subsequent 
dust grains have to be taken into account. If we denote the 
total number of scattering events with M, the velocities of 
the dust grains with and the propagation directions with 
ki (for i — 1 . . . M), we obtain 



■ feo + v d% ■ 



(10) 



In general, therefore, the inclusion of kinematical informa- 
tion into radiative transfer problems is very complex. For our 
problem of dust attenuation in elliptical galaxies, however, 
we can make one assumption that reduces the complexity 
considerably: the summation in formula (|lo|), i.e. the com- 
mon contribution of the dust velocities, is in general neg- 
ligible with respect to the first term, the contribution of 



the stellar velocity. Several arguments support this assump- 
tion. [1] It is reasonable to assume that the dust grains have 
smaller velocities than the stars in elliptical galaxies. In- 
deed, if the dust grains were to have large velocities, they 
would collide and heat up. However, the lion's share of the 
dust in elliptical galaxies is assumed to be cold (Fich & 
Hodge 1993; Goudfrooij & de Jong 1995; Wiklind & Henkel 
1995; Merluzzi 1998). [2] Scattering off interstellar dust is 
generally anisotropic, with a larger probability for forward 
scattering (see Section 1.2.2] ). In Appendix B we show that 
the anisotropic nature of scattering contributes to reducing 
the importance of the dust grain velocity terms in equation 
(0). [3] There is no reason why the velocity of dust grains 
in an elliptical galaxy would have a preferential direction. 
The individual terms in the summation in formula (|l^) will 
therefore have a random sign, and the summation will in the 
mean be washed out by multiple scattering. 

As a consequence, we can write equation ( |lo| simply as 
Vii — v* ■ feo. Each photon hence carries the velocity com- 
ponent of the star in the direction of the emission of the 
photon. 



3 DESCRIPTION OF THE MONTE CARLO 
ROUTINE 

3.1 Basic characteristics 

We constructed a Monte Carlo code to solve our radiative 
transfer problem. Usually, the main argument against Monte 
Carlo methods is that it is computationally rather expen- 
sive compared to other methods (see e.g. Baes & Dejonghe 
2001a). It has other advantages, however, which make it very 
competitive, in particular in an era when CPU time is not 
the most stringent limitation anymore. Important qualities 
of the Monte Carlo method are the possibility for a proper 
error analysis and a very wide flexibility. Nice examples of 
this flexibility are the ability to handle arbitrary geome- 
tries (Wolf, Fischer & Pfau 1998; Wood & Reynolds 1999; 
Gordon et al. 2001), the polarization of the scattered radia- 
tion (Code & Whitney 1995; Bianchi, Ferrara & Giovanardi 
1996), the dumpiness of the interstellar medium (Witt & 
Gordon 1996; Bianchi et al. 2000), and the self-consistent 
heating and re-emission of the absorbed radiation (Wolf, 
Henning & Stecklum 1999; Bianchi, Davies & Alton 2000; 
Misselt et al. 2001). We make use of another possibility of 
the Monte Carlo method: its ability to include kinematical 
information in an elegant way. This possibility has recently 
been explored by Matthews & Wood (2001), who studied the 
effect of dust attenuation on the observed rotation curve in 
the spiral galaxies. Our modelling opens up the possibility 
to construct all LOSVDs and hence to investigate the entire 
observed kinematical structure of galaxies. 

The principles of the Monte Carlo technique are out- 
lined in detail by various authors (Cashwell & Everett 1959; 
Mattila 1970; Witt 1977; Yusef-Zadeh, Morris & White 1984; 
Fischer, Henning & Yorke 1994; Bianchi et al. 1996). Basi- 
cally, the routine consists of following the individual trajec- 
tory of a very large number of photons through the galaxy. 
A photon's history is given by a number of quantities such 
as the position and propagation direction at birth, the dis- 
tance the photon travels before it interacts with a dust grain, 
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the kind of this interaction, etc. Each of these quantities 
is described statistically by a random variable, taken from 
a particular probability density. Our approach is based on 
Witt (1977) and Bianchi et al. (1996). It includes the use 
of a continuous three-dimensional cartesian reference sys- 
tem (no grid). Furthermore, we use the classical tricks to 
optimize the routine: we assign a weight to each photon in 
order to avoid the loss of photons due to absorption, and 
we apply the forced first scattering mechanism to improve 
the statistics of the scattered radiation (Cashwell & Everett 
1959; Witt 1977). The novelty of our Monte Carlo code is 
that it can calculate both photometric and kinematic data: 
not only the light profile, but also the LOSVDs and the 
projected kinematic moments. It simultaneously calculates 
these data in three modes: without dust attenuation, with 
only absorption taken into account and with dust attenua- 
tion fully taken into account. It is a monochromatic code, 
i.e. it calculates the observables in one single wavelength. 
If desired, the wavelength dependence of the results can be 
investigated by running the code various times at different 
wavelengths. We implemented the central wavelengths of the 
optical U, B, V, R, I and the near-infrared J, H and K 
broadbands. 

3.2 Calculation of the light profile 

The light profile of the galaxy is constructed by simulating 
an imaging process. If a photon leaves the galaxy, we record 
the position of the last scattering event and the final prop- 
agation direction of the photon. Because in spherical sym- 
metry, a path can be determined by its projected radius, 
we only need to record the projected radius of the photon's 
final path. The registration of the photons is done by con- 
structing a histogram of the photons leaving the galaxy as 
a function of 7?, where the contribution of each photon is 
measured by its weight (see Witt 1977). 



conditions, we follow the photon until the exit conditions 
are satisfied, and calculate the projected radius of the final 
line of sight. The photon will then be stacked, according to 
its weight, in the appropriate bin in a two-dimensional ar- 
ray with axes R and vn . This process is repeated this for 
a large number of photons and two-dimensional histograms 
are formed. 

This method can be optimized in two ways. Firstly, we 
see that the velocity vector v* does not show up in the Monte 
Carlo procedure itself, but only in the classification process. 
Moreover, we actually don't need the entire information con- 
tained in the velocity i>»: we only need the component v\\ of 
the velocity in the direction feo. Therefore it is sufficient to 
generate a random vn from the spatial LOSVD <fi(ro, feo, vu ), 
which represents the probability density of line-of-sight ve- 
locities for a star at position ro in the direction feo. Instead of 
generating a three-dimensional velocity vector for each pho- 
ton in the beginning of the Monte Carlo cycle, it is hence 
sufficient to generate one single line-of-sight velocity com- 
ponent. This saves many random number generations and 
distribution function evaluations, which are costly processes. 

A second optimization eliminates the generation a vn 
altogether. Instead of generating one single random v« , we 
can, for each bin in the velocity direction, calculate the prob- 
ability that the line-of-sight velocity will fall in that bin. If 
the boundaries of a bin are given by and v\\j, this 

probability equals 

/""Id 

/ (j>{ri,k ,v\i)dv n 
Jv \\,i-i 

^ 1 <ft(ri,fc ,a| M -i) + <f)(r 1 ,k ,v ]lj ) 
~ 2 v \\,j- v hj-i 

For each photon, we assign a weighted value to each velocity 
bin corresponding to the final projected radius, which gives 
better statistics than dropping the photon in a single bin. 



3.3 Calculation of the LOSVDs 

The Monte Carlo method provides an elegant possibility to 
include kinematical information. Indeed, we argued that the 
velocity information that a photon carries is essentially the 
velocity component of the star in the direction of the emis- 
sion. This extra information can easily be included in the 
Monte Carlo routine. The only extension is that we do not 
only have to generate an initial position and propagation di- 
rection for each photon, but also a stellar velocity. The initial 
position and the stellar velocity have to be extracted in ac- 
cordance with the phase space distribution function F(r, v), 
which represents the probability of finding a star at the posi- 
tion r with velocity v. A practical way to do so, is following 
the strategy of Wybo & Dejonghe (1996), who generated 
sets of phase space coordinates of stars to create iV-body 
representations of globular clusters. First, a random posi- 
tion ro is determined from the stellar emissivity, and then 
a random stellar velocity v, is generated from the three- 
dimensional probability density p(v ) — C F(ro, v )/l(ro), us- 
ing the acceptance-rejection method (Abramowitz & Stegun 
1965; Press et al. 1989). Next, we determine a random emis- 
sion direction feo, and determine the velocity component of 
v, in the direction feo, i.e. v\\ = v, ■ feo. Given these initial 



3.4 Calculation of the projected kinematics 

Of course, we are not able anymore to write down a direct 
connection between the intrinsic moments of the distribution 
function and the projected moments on the plane of the sky, 
contrary to the absorption-only case (Paper I). Instead, we 
will calculate the projected moments directly from the ob- 
tained LOSVDs. The projected kinematics we consider are 
the mean projected velocity v p , the projected velocity dis- 
persion (j p , the skewness ks, the kurtosis fc.4, and the lowest 
order Gauss- Hermite coefficients /13 and /14. Note, however, 
that for spherically symmetric non-rotating galaxies all odd 
moments vanish, i.e. v p = £3 = /13 = 0. 

3.5 Error analysis 

Obviously, a reliable error analysis is necessary to estimate 
the accuracy of the results. But it will be even more impor- 
tant in a later stage. Indeed, we are presently incorporating 
this radiative transfer mechanism into the QP program (De- 
jonghe 1989), a code for the dynamical modelling of gravitat- 
ing systems. Therefore, it is necessary that we can dynami- 
cally evaluate the error bars on any observable: the program 
has to be able to calculate any observable with an preset 
accuracy, required by either the user or another part of the 
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code itself. Monte Carlo methods satisfy this need: they can 
keep on adding photons until any accuracy requirement is 
satisfied. 

The determination of the errors is easy for observables 
that are directly proportional to the number of photons de- 
tected in the bins, e.g. the intensity. Because of the Poisson 
character of the noise, the error bars on these quantities can 
be directly calculated from the square root of the number of 
photons. The observed kinematics we calculate, however, are 
not directly proportional to number of photons obtained in 
the bins, such that we cannot adopt this simple procedure. 
Therefore, we estimate the error on our observables using 
an alternative method. Assume we have calculated a value 
/ for a given observable (either photometrical or kinemati- 
cal) after a run with TV photons. We divide the total number 
of photons in M subsets, and for each of these subsets, we 
calculate the corresponding observable fi. The uncertainty 
on / is estimated as, 



1 



1 1/2 



fY 



(12) 



Typically, we set M — 100, and we adopt a minimum num- 
ber of 10 6 photons, i.e. 10 4 photons per subset. 



4 THE GALAXY MODELS 

To investigate the effects of dust attenuation on the light 
profile and the observed kinematics, we will adopt a set of 
dusty galaxy models, similar to those in Paper 1 and Pa- 
per II. They are spherically symmetric, and consist of a 
stellar and a dust component. 



4.1 The stellar component 

For the stellar component we use a self-consistent Plum- 
mer model (Plummer 1911; Dejonghe 1987), defined by the 
potential-density pair 

GM 



p{r) 



\/r 2 + c 2 



3M 
4tv 



1 + - 



-5/2 



(13a) 



(13b) 



We adopt a core radius c = 5kpc, a total mass M — 6 x 
10 11 Mq and a (constant) mass-to-light ratio T = 4 Yq, such 
that the total emitted luminosity equals C = 1.5 x 10 11 L©.^| 
As in Paper I, we consider different Plummer models, each 
with a different value for the parameter q, which describes 
the internal anisotropy of the stellar orbits. We adopt three 
different models: a radial model (q — 2) , the isotropic model 
(q = 0) and a tangential model (q — —6). These are the 
same values as in Paper I, except that we replace the radial 
model with q = 1 by the model with q = 2. The reason 
for preferring q — 2 above q = 1 has a computationally 
nature: the spatial LOSVDs can be calculated analytically 

3 These parameters differ from those adopted in Paper II, which 
contains an error: the adopted mass is inconsistent with the ve- 
locity dispersion. When we refer to results of Paper II, these are 
scaled to the parameters adopted in this paper. 



for q = 2, whereas this is not possible for the q 
(Appendix A). 

4.2 The dust component 



1 model 



As shown in Section 2.1, a complete characterization of the 



dust component requires a specification of the opacity n(r), 
the albedo ui(r) and the scattering phase function $(r, k, fc'). 
The opacity sets the total amount of dust and its spatial 
distribution, whereas the other two quantities describe the 
optical properties at a given distance to the galaxy centre. 
In order to limit the number of parameters in our model, 
we will assume that the optical properties of the dust grains 
are the same all over the galaxy, such that the albedo and 
the phase function become independent of r. We will return 
to this assumption in Section 5.3.3. 

4-2.1 The dust distribution 

We adopt the same family of dust components considered 
in Paper I and Paper II, characterized by the dust number 
density 

■x/2 



n(r) = no I 1 + 



(14) 



where no is the dust density at the center of the galaxy, and 
c and a determine the actual shape of the dust distribution. 
Because the effect of varying a and c are roughly similar, we 
restrict ourselves to varying the dust exponent a, and we fix 
the dust core radius to c = 5 kpc, equal to the stellar core 
radius. The variation of a still allows us to explore a large 
variety of star-dust geometries: small values of a correspond 
to centrally concentrated distributions, whereas for large a 
the distribution is shallower (Paper I). 

At a fixed wavelength A, the opacity function n\(r) can 
be derived from the dust density by multiplication with the 
(wavelength-dependent) dust cross-section s\. To character- 
ize the dust distribution at a given wavelength, we thus need 
three parameters: the dust exponent a, the central dust 
number density, and the dust cross-section s\. It is more 
convenient, however, to use another, equivalent, set of pa- 
rameters. Indeed, if the opacity function is known at one 
particular wavelength, it can be determined at any other 
wavelength. Usually, the F-band is used as the reference 
wavelength. We can then write «a (r) ~ X\ ny (r) , where 
X\ = sx/sv is the extinction coefficient, relative to the V 
band. It is then customary to introduce the total F-band 
optical depth as[|] 



TV = 



Ky(r) dr. 



(15) 



4 In this definition, ry is defined as the optical depth from the 
center of the galaxy to the edge. We want to stress that this 
definition differs a factor 2 with the definition of the optical depth 
we adopted in Paper I and Paper II, where we defined ry as the 
total optical depth along the central line of sight, i.e. from one 
edge to the other. The definition of ry we adopt here, however, 
is more commonly used, and it will allow the reader to link the 
obtained results more easily to those found in the literature (e.g. 
Witt et al. 1992; Wise & Silva 1996). When we refer to the results 
of Paper I and Paper II, we will adopt the convention used here. 
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Table 1. The adopted optical properties of the dust grains, at 
the central wavelengths of the standard optical and near-infrared 
wavebands. Tabulated are the relative extinction coefficient X\, 
the scattering albedo lux and the asymmetry parameter g x . 



band 


A (nm) 


x x 


MX 


9a 


U 


360 


1.52 


0.63 


0.65 


13 


440 


1.32 


0.61 


0.63 


V 


550 


1.00 


0.59 


0.61 


R 


700 


0.76 


0.57 


0.57 


1 


850 


0.48 


0.55 


0.53 


J 


1250 


0.28 


0.53 


0.47 


11 


1650 


0.167 


0.51 


0.45 


K 


2200 


0.095 


0.50 


0.43 



If we use the triplet (tv , a, X\) as parameters, the opacity 
at an arbitrary wavelength A can be written as 

/2 

(16) 



9 r - 



X\TV 



x/¥r(^i) c 



T 

1 + - 

& 



4.2.2 The phase function 

We can safely assume that the phase function does not de- 
pend independently on the four variables (k,k'), but only 
on the angle between the two directions k and k! . Set- 
ting k ■ k' = cos a, we can write the phase function as 
$(cosa). Probably the most widely adopted phase function 
that can describe anisotropic scattering is the one named af- 
ter Henyey & Greenstein (1941). It has been derived empir- 
ically, as a description of the scattering of light in reflection 
nebulae in the Galaxy, and has a simple analytical form, 



$(cosa) = 



(l+g 2 -2g cos a) 3 / 2 ' 



(17) 



This family of phase functions contains the so-called asym- 
metry parameter g, which is a measure for the anisotropy 
of the scattering. In particular, for g — the scattering is 
isotropic, whereas for g = 1 the scattering is completely 
forward. 



4-2.3 The optical properties 

From the previous subsections, we see that, in total, we need 
five parameters to describe a dust component at a given 
wavelength A: the V-band optical depth tv, the dust expo- 
nent a, the relative extinction coefficient X\, the scattering 
albedo lu\ and the asymmetry parameter g\. In principle, 
all of these parameters could be considered as free parame- 
ters. In order to limit the degrees of freedom, however, we 
decided not to consider the optical properties of the dust 
as a set of free parameters, but to use a fixed set of dust 
parameters. The values we use are listed in Table [l]. They 
are taken from Gordon, Calzetti & Witt (1997), who derived 
them from an interpolation between a large set of empirical 
data of typical Milky Way dust. Notice in particular that 
the scattering albedo does not vary greatly within the op- 
tical and near-infrared regime. It is always slightly greater 
than 0.5, i.e. the probability for absorption is always slightly 
smaller than the probability for scattering. 



4.3 A template model 

Our models contain three parameters, the orbital structure 
parameter q, the V-band optical depth tv and the dust ex- 
ponent a. By varying these, we can describe a set of dusty 
elliptical galaxy models with a large variety in internal dy- 
namics, dust content and star-dust geometry. We apply the 
same strategy as in Paper I to present the results of our 
modelling. First, we keep the geometry of the dust compo- 
nent fixed, and we investigate how dust distributions with 
various optical depths affect the light profile and observed 
kinematics of our three Plummer models. Next, we keep the 
optical depth fixed and consider a set of star-dust geometries 
by varying the dust exponent a. 

We therefore require a template model, where we can 
vary one parameter while keeping the other one fixed. For 
the optical depth of this template model, we choose tv = 1, 
a value that could be in agreement with both the far-infrared 
emission and the colour gradients in elliptical galaxies (Witt 
et al. 1992; Goudfrooij & de Jong 1995; Wise & Silva 1996). 
For the dust geometry of the template model, we adopt a 
modified Hubble profile, characterized by setting a = 3. It 
is shallower than the stellar distribution, which seems indi- 
cated by colour gradient models (Wise & Silva 1996). As a 
note we would like to stress that both the optical depth and 
geometry of the diffuse dust component in elliptical galaxies 
are very poorly known, apart from the solid result that FIR 
emission indicates that it is present. We return to this issue 



in Section 6.2. 

Finally, we need to choose a template wavelength to 
present our results. Because the present paper focuses on the 
effects of dust attenuation on the observed kinematics, and 
because most of the kinematical observations are presently 
conducted at optical wavelengths, we will work in the V 
band, unless mentioned otherwise. 



5 RESULTS 

5.1 Putting the method to the test 

The Monte Carlo routine yields results in three modes, 
the nodu, abso and dust modes. Whereas we are of course 
mainly interested in the dust results, the other data sets are 
useful as a check for the accuracy of our results. If dust is not 
taken into account, the light profile and the projected kine- 
matics can be calculated analytically for the galaxy mod- 
els we consider (Dejonghe 1987). If only dust absorption is 
taken into account, they can be calculated through one sin- 
gle quadrature (Paper I). In Fig. |^ we compare the results 
of our Monte Carlo code with the corresponding analytical 
results, for two different values of N, the total number of 
emitted photons. Even for N as low as 10 6 , the minimum 
number of photons we consider, the analytical results are 
very well reproduced, and are everywhere within the error 
bars. 



5.2 The light profile 

5.2.1 Dependence on the optical depth 

In Fig. |^ we demonstrate how the light profile is affected 
when both absorption and scattering are included. In the 
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Figure 2. A comparison of the results of the Monte Carlo code in nodu mode (grey curves with error bars) and the corresponding 
analytical results (black lines). The upper and lower panels correspond to different values for the total number of emitted photons. 
Shown are the surface brightness profile, the projected velocity dispersion and the /14 profile, for the three different orbital modes 
considered. 



center of the galaxy (where the lines of sight contain most of 
the dust), the attenuation is strongest, and the attenuation 
decreases as one goes to the outer regions. As a consequence, 
also the apparent luminosity decreases, and the apparent 
size of the core (as measured by the effective radius R e ) 
increases as a function of Ty. 

Within the first few kpc, dust attenuation has at first 
order the same effect on the light profile as absorption alone. 
Roughly, the effects of a attenuating dust component with 
optical depth tv can be approximated by a purely absorb- 
ing dust component with effective optical depth tv/2. This 
approximation has a natural explanation: it follows from as- 
suming that the scattering is purely forward. Indeed, the 
phase function corresponding to forward scattering is a sim- 
ple Dirac delta function, <I>(fc,fc') = An 8(k — k'), and sub- 
stituting this into the RTE (g), we find 

^-{r, k) = e(r) - (1 - w) K(r) I(r, k). (18) 
as 

This radiative transfer problem is completely analogous to 
a radiative transfer problem where only absorption is ac- 
counted for, but where the optical depth of the dust com- 
ponent is diminished with a factor 1 — w w 5. This ap- 



proximation is an appealing way to estimate the effects of 
scattering without elaborate and costly radiative transfer 
calculations. An argument in favour of this approximation 
is that multiple scattering events tend to wash out the effects 
of anisotropy of the scattering phase function. Hence, in me- 
dia with a large opacity, any phase function can in principle 
be adopted, including the degenerate one corresponding to 
forward scattering (e.g. Di Bartolomeo, Barbaro & Perinotto 
1995). 

However, the effects of scattering should never be under- 
estimated. In a previous study on the RTE in plane-parallel 
geometry, we investigated several ways often used in the lit- 
erature to approximate scattering (Baes & Dejonghe 2001b). 
One of our main conclusions was that none of these methods 
provides a satisfactory approximation to the exact solution. 
The reason is that the physical process of scattering has 
a completely different nature than absorption, because it 
changes the path along which the photons propagate. The 
effect of scattering is that photons have a preferential direc- 
tion in which to leave the galaxy. In a plane-parallel geome- 
try, photons prefer the face-on direction to leave the galaxy, 
because when a photon is scattered into that direction, its 
chances to leave the galaxy are larger (the optical depth 
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Figure 3. The effects of dust attenuation on the light profile, as a function of the optical depth. The left panel shows the attenuation 
profile of the galaxy, i.e. the fraction of the light output attenuated by the dust in magnitude units. It is shown for absorption only 
(dashed lines), and for full attenuation (solid lines), for optical depths ry = 0, ^, 1 and 2 - only the extreme ones are labelled. The middle 
plot shows the cumulative luminosity function, i.e. the fraction of the total emitted luminosity detected inside an aperture of radius R. 
The plots on the right-hand side show the effect of dust attenuation on the total attenuation and the effective radius, as a function of 
optical depth. 



along the path is shorter) than when it is scattered into 
inclined directions. An analogous reasoning holds for spher- 
ical geometry, where photons will generally prefer large pro- 
jected radii to leave the galaxy.^ As a consequence, the net 
effect of scattering is that photons are scattered out of lines 
of sight close to the center of the galaxy, into lines of sight 
with a large projected radius. In the outer regions of the 
galaxy, these extra photons will reduce the loss of radiation 
due to absorption. At very large projected radii, the attenu- 
ation will even be negative, i.e. the galaxy will even appear 
brighter than when dust extinction is not taken into account 
(Fig. ^|). These results are in agreement with those found by 
Wise & Silva (1996), and they will be very important when 
we investigate the full effects of attenuation on the observed 
kinematics. 

5.2.2 Dependence on the dust geometry 

In Fig. ^]we show the effects of varying the dust exponent on 
the light profile of our galaxy models. In the central regions 

5 At any position in a spherical galaxy, the optical depth is obvi- 
ously smallest along the path in the radial direction, i.e. the path 
directed away from the center of the galaxy. One would there- 
fore be inclined to think that the photons will prefer to leave the 
galaxy through the central lines of sight. However, the connection 
between lines of sight and directions is not one-to-one: to every 
possible line of sight a photon can be scattered into (i.e. to ev- 
ery R ^ r), correspond two directions, one directed towards the 
center of the galaxy and one towards the edge of the galaxy. The 
probability that a photon, scattered at r into a line-of-sight R, 
will leave the galaxy, equals the weight function K(R, r), defined 
in Section 2 of Paper I. For fixed values of r, the function IC(R, r) 
is an increasing function of R, hence having its maximal value 
at R = r. Photons will hence on average more easily leave the 
galaxy at large projected radii. 



of the galaxy, the combined effects of scattering and absorp- 
tion on the light profile can roughly be approximated by pure 
absorption with an effective optical depth r{> ft = (1 — tv) Ty. 
Therefore, the dependence on the dust exponent are sim- 
ilar to those when absorption only is taken in to account 
(Paper I): the attenuation is stronger for extended dust dis- 
tributions (a small) than for centrally concentrated ones (a 
large). Because the central regions emit most of the light 
and therefore dominate the total observed luminosity, also 
the total attenuation will decrease as a function of a. 

At large projected radii, the dust has another effect on 
the light profile: the net effect is that photons scattered into 
these lines of sight cause a negative attenuation, i.e. the 
galaxy appears brighter. This effect should be stronger for 
extended dust distributions, because these distributions im- 
ply more dust in the outer regions of the galaxy, and hence 
an enhanced probability to be scattered into the outer lines 
of sight. This can indeed be observed in Fig. |^. In partic- 
ular, this negative attenuation is nearly non-existent for a 
model where dust and stars have the same spatial distri- 
bution (a — 5), whereas it is clearly noticeable if the dust 
distribution is shallower than the stellar distribution. 

5.3 The observed kinematics 

5.3.1 Dependence on the optical depth 

In an optically thin galaxy, the LOSVD is formed by sum- 
ming the contribution of the line of sight velocities of all 
stars that are situated along that line of sight. If absorption 
is included in the projection process, still the same stars on 
the line of sight contribute to the LOSVD, but the contribu- 
tion of each star is weighted by the amount of starlight that 
is able to survive the absorption and reach the observer. The 
net effect is that, for a given line of sight, the stars at the 
outer parts of the line of sight contribute relatively more to 
the LOSVD than the stars in the central parts (Paper I). 
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Figure 4. The effects of dust attenuation on the light profile, as a function of the dust geometry. The four panels in this figure are 
analogous to those in Fig. ^. The models shown here all have Ty = 1, but different values for a. The values shown are a = 1.25, 2, 3 and 
5. The black dotted lines in the left and middle panels correspond to the model without dust attenuation. The grey dotted lines on the 
right-hand side represent models with only absorption taken into account, as in Fig. 0. 



Because the largest line-of-sight velocities along a given line 
of sight are usually found around the tangent point, the 
LOSVDs will be biased towards smaller velocities. In par- 
ticular, the projected velocity dispersion generally decreases 
if absorption is taken into account. Only for the outer lines 
of sight of galaxies with a very radially anisotropic orbital 
structure, the largest line-of-sight velocities are found in the 
outer regions. Indeed, the stars at the tangent point have 
small line-of-sight velocities, because their (radial) orbits are 
nearly perpendicular to the line of sight. We showed in Pa- 
per I that the effects of dust absorption on the LOSVDs are 
only considerable for large optical depths (tv > 5), which 
are probably not appropriate for elliptical galaxies. 

If scattering is taken into account, the situation changes 
drastically. The net effect of scattering, at least for the light 
profile, is that photons are taken away from the central lines 
of sight and sent into lines of sight at larger projected radii. 
Now consider a photon that is emitted by a star near the 
center of the galaxy, and that, after one or several scattering 
events, propagates towards the observer at a large projected 
radius. Although this photon will contribute to the LOSVD 
at this large projected radius, it carries the velocity informa- 
tion from the emitting star. Notice that this star does not 
physically belong to that line of sight. Hence, when scat- 
tering is taken into account, the LOSVDs aren't LOSVDs 
anymore in the strict meaning of the word: the LOSVD at 
a certain line of sight can contain information of stars at 
totally different lines of sight. More generally, every single 
star in the galaxy will contribute to every single LOSVD, 
whereby its contribution will be weighed by the number of 
photons that leave the galaxy along that line of sight. 

Beside the photons that disappear from the line of 
sight due to attenuation, we hence also have to account 
for the photons scattered into the line of sight, which 
contribute the additional kinematical information from 
stars that physically don't belong there. How this process 
affects the LOSVDs is more complex than the effects due to 
the photons taken away from the line of sight. In particular, 



this effect is different for lines of sight that pass through 
the center of the galaxy (R « 0) and lines of sight at larger 
projected radii (R c). 

If scattering is not taken into account, the LOSVD at 
the central line of sight R — only contains information 
on the radial velocity component of stars, because along the 
central line of sight Vu = v r . Because the vast majority of the 
stars along this line of sight reside in the central regions, the 
central LOSVDs will be dominated by the distribution of ra- 
dial velocities in the galaxy center. If scattering is included, 
the central line of sight will also contain photons scattered 
into it, that would normally leave the galaxy at a larger pro- 
jected radius. The vast majority of these photons will also 
originate from near the central region, but the line-of-sight 
velocities carried by them are not necessarily the radial ve- 
locity components of the stars that emitted them. These 
photons contaminate the LOSVD with tangential velocities: 
instead of containing pure radial velocity information, it will 
reflect a mix of radial and tangential components. As a con- 
sequence, the effect of scattering will depend on the ratio of 
radial to tangential velocity components in the central re- 
gions of the galaxy, i.e. on the orbital structure. In the panels 
on the middle row in Fig. ^, we plot the effect of increasing 
optical depth on the central projected velocity dispersion. 

For isotropic galaxies, radial and tangential velocities 
are in balance throughout the galaxy, such that the photons 
scattered into the central lines of sight will hardly affect 
the LOSVD. The total effect of attenuation on the central 
LOSVDs will be dominated by the absorption effect, i.e. a 
bias towards smaller line-of-sight velocities. In particular, 
the central projected dispersion decreases if dust is taken 
into account, in a very similar way as when only absorption 
is taken into account. 

In a radially anisotropic galaxy, stars have on average 
a larger radial than tangential velocity component. Stars 
scattered into the central line of sight that contribute part 
of their tangential velocity, will therefore bias the LOSVD 




Figure 5. The effects of dust attenuation on the observed velocity dispersion, as a function of the optical depth. In the top panels, the 
projected velocity dispersion profiles are plotted, for different values of Ty. We have used the same values of ry and the same layout as in 
Fig. ^. The left, middle and right columns correspond to our radial, isotropic and tangential models respectively. The three panels on the 
middle row show the central projected velocity dispersion as a function of Ty. The values are normalized by cr P: o, the central projected 
velocity dispersion in an optically thin galaxy. The solid lines represent the effect of dust attenuation, whereas in dashed lines we show 
the effects of absorption only. Analogously, the three bottom panels show the effects of dust attenuation, respectively absorption, on the 
normalized projected velocity dispersion at R = 25kpc. 



towards smaller velocities. Because the effect of absorption 
is also a bias towards smaller velocities, the total effect of 
attenuation on the central LOSVDs of radially anisotropic 
galaxies is to turn them more peaked, i.e. to decrease the 
central projected dispersion. The strength of the effect is 
stronger than the effect of absorption alone. 

On the contrary, for tangentially anisotropic galaxies, 
the tangential velocity of the stars generally exceeds their 
radial velocity. As a consequence, if the central LOSVD 
is contaminated with tangential velocity information, it 
will be biased towards larger line-of-sight velocities. Two 
processes hence affect the LOSVD in an opposite way: 
absorption (and scattering of photons out of the line 
of sight) favors smaller velocities, whereas scattering of 
photons into the line of sight biases the LOSVD towards 
larger velocities. The way the central LOSVDs are affected 
depends on which of both mechanisms is stronger. For the 



modest optical depths appropriate in elliptical galaxies, 
the scattering effect dominates, and the central projected 
dispersion increases due to scattering. 

For lines of sight at large projected radii, the situation 
is completely different. We argued that the net effect of scat- 
tering is that photons are scattered from the central lines of 
sight into the outer lines of sight. These scattered photons 
will strongly affect the observed kinematics. Indeed, the ma- 
jority of these scattered photons are emitted by stars in the 
central regions of the galaxy. They carry along the kinemat- 
ical information of the stars that emitted them, i.e. the typi- 
cal velocities appropriate in the central regions of the galaxy. 
These are on average larger than the typical line-of-sight ve- 
locities in the outer regions of the galaxy, where the kinetic 
energy of the stars is much lower. As a consequence, they will 
contaminate the LOSVDs with high velocities. In particular, 
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Figure 6. The effects of dust attenuation on the observed LOSVD shape, parametrized by the kurtosis fej and Gauss-Hcrmite /14 
coefficient, as a function of the optical depth. As in Fig. ^| the left, middle and right columns correspond to the radial, isotropic and 
tangential models respectively. The adopted optical depths and layout are similar to those used in Fig. H. 



it should be noted that lines-of-sight velocities will be ob- 
served in the LOSVDs, that would normally be impossible 
at these projected radii. If scattering is not taken into ac- 
count, the maximal line-of-sight velocity of the LOSVD at a 
projected radius R is the escape velocity v csc (R) — ^2tp(R), 
where ip(R) is the potential at the tangent point r = R. If 
scattering is taken into account, the photons observed at the 
line-of-sight R can originate from stars in the center of the 
galaxy where the line-of-sight velocity can be larger that the 
escape velocity at 7?. These large line-of-sight velocities will 
cause "forbidden" high- velocity wings in the outer LOSVDs. 
As a result, the projected dispersion at large projected radii 
will increase. 

In the top panels of Fig. ^, we plot the projected veloc- 
ity dispersion profiles for various values of the optical depth. 
Clearly, the velocity dispersion increases significantly in the 
outer regions. The strength of this increase depends on the 
orbital structure of the galaxy. Radially anisotropic models 
are more strongly affected than tangential ones, because 
their (optically thin) outer LOSVDs are more strongly 
peaked, and thus more vulnerable to the contribution of 
photons from high-velocity stars in the center. This is 
clearly shown in the bottom row of Fig. H, where the effect 



of dust attenuation on the projected velocity dispersion at 
large projected radii is shown as a function of the optical 
depth. The effects are significant: for example, for an optical 
depth of unity, the projected dispersion increases with more 
than 40 per cent for the radial model. 

Of course, the LOSVDs are not completely determined 
by the projected velocity dispersion alone, which just gives 
a measure for the broadness of the LOSVD. The extra in- 
formation contained in the (symmetric) LOSVDs, the ac- 
tual shape of the LOSVDs, can be represented by either 
the kurtosis or the /14 parameter. In the top panels of Fig. 
^j] we show how dust attenuation affects the shape of the 
LOSVDs of our Plummer galaxies, as quantified by the kur- 
tosis. In the inner regions, the effect is fairly small. As the 
effect on the projected dispersion, its sign is dependent on 
the orbital structure: for the radial and isotropic models, the 
kurtosis increases with optical depth, whereas it decreases 
for the tangential model. At large projected radii, the kurto- 
sis increases spectacularly if dust attenuation is taken into 
account, and again, the effect is much stronger for radial 
than for tangential models. The reason for this increase are 
the high- velocity wings in the LOSVDs. 
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Figure 7. The effects of dust attenuation on the observed kinematics, as a function of the dust geometry. Shown are the projected 
dispersion profile and the /14 profile, for different values of the dust exponent a. The presented models and the layout are the same as 
in Fig. H. As in the previous figures, the three column correspond to the three different Plummer models. 



The bottom panels in Fig. ^| show the effects of dust 
attenuation on the shape of the LOSVDs of our models, as 
quantified by the /14 parameter. The effects are comparable 
to those on the kurtosis. This comes as no surprise, because 
kurtosis and /14 are proportional to first order (van der Marel 
& Franx 1993). 

5.3.2 Dependence on the dust geometry 



In Section 5.2.2 , we showed that the contribution of scat- 
tered photons to the surface brightness at large radii de- 
pends rather critically on the dust geometry. If, on the one 
hand, the dust distribution is shallower than the stellar dis- 
tribution, the contribution of these photons is important, 
and they even cause a negative attenuation at the outer lines 
of sight. If, on the other hand, the dust follows the stellar 
distribution, the effects of scattering are nearly negligible. 
Because we showed that these scattered photons affect the 
observed kinematics rather strongly, it can be expected that 
the effect of dust attenuation on the observed kinematics 
will also be very sensitive to the dust geometry. 

In Fig. Q we illustrate how the observed kinematics de- 
pend on the dust geometry. We find indeed that the kine- 



matics are much more affected for extended dust distribution 
than for centrally concentrated ones. In particular, when the 
dust has the same spatial distribution as the stars, the ef- 
fects of scattering on the observed kinematics are nearly 
negligible. On the contrary, when the dust density decreases 
very slowly, the outer regions have a large dust-to-stars ra- 
tio, such that the photons scattered into these lines of sight 
form a large fraction of the total number of photons that 
contribute to the LOSVDs. The high-velocity stars gradu- 
ally contribute more as the dust geometry becomes more ex- 
tended. As a consequence, the projected velocity dispersion 
profile increases significantly with decreasing a at large pro- 
jected radii, and also the effect on the /14 shape parameter 
depends strongly on the dust geometry. The dust distribu- 
tion is clearly an important parameter in our models. 

5.3.3 The influence of optical property gradients 

From the previously obtained results, we know that the ef- 
fects of dust attenuation on the observed kinematics are due 
to photons emitted by high-velocity stars in the centre of 
the galaxy, scattered in the outer regions on lines of sight at 
large projected radii. The strength of these effects will hence 
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Figure 8. The effect of a dust albedo gradient on the attenuation curve and the observed kinematics, more precisely the projected 
dispersion and /14 profiles. The adopted model is the isotropic template model in the V band, but with a variable dust albedo of the 



form (J19p , with luq = u)y = 0.59, 
model without dust attenuation. 



and three different values of LUoo, being 0.79, 0.59 and 0.39. The black dotted curves correspond to the 



depend on the probability that a photon from a high veloc- 
ity star will be scattered onto outer lines of sight. There are 
various factors that can contribute to this probability. We 
already encountered two of them: the larger the dust con- 
tent of the galaxy (high r) and/or the more extended the 
dust distributions (small a), the greater the probability for 
scattering at large radii, and hence the stronger the effects 
on the observed kinematics. But also the optical proper- 
ties of the dust can contribute to the number of scattering 
events that send photons from the inner lines of sight to 
outer lines of sight. Imagine, for example, that the scatter- 
ing would be completely forward at large radii, or that the 
scattering albedo would become negligible in the outskirts 
of the galaxy. In either case, the probability that a photon 
emitted in the centre would reach the observer along a line 
of sight at large projected radii would be very small, such 
that the observed kinematics would hardly be affected by 
dust attenuation. 

So far, we have adopted the assumption that the optical 
properties of the dust are equal all over the galaxy. However, 
this assumption might not always be satisfied in real ellip- 
tical galaxies. Indeed, gradients in e.g. the metallicity, the 
stellar radiation field or the X-ray luminosity density of a 
galaxy can cause systematic changes in the size distribu- 
tion and/or chemical composition of the dust grains, which 
can result in gradients in their optical properties. In our own 
Galaxy, the physical properties of interstellar dust have been 
found to vary substantially in different environments (Witt, 
Bohlin & Stecher 1984; Mathis & Cardelli 1992). 

With our Monte Carlo technique, we can easily include 
optical properties gradients in our models. The strongest 
effect can be expected for a gradient in the scattering albedo, 
because w(r) is directly related to the number of scattering 
events at a distance r. We add to the template model of 
Section 4.3 an albedo gradient of the form 



lo c 2 + io x r 2 



(19) 



c 2 _|_ r 2 ' 

which changes smoothly from u>o in the centre to luoo at 



large radii. In Fig. [8] we show the effect of such a gradient 
on the a p and fu profiles, with the central albedo wo the 
value from Tab. [I], and a number of different values for lo x . 
This figure demonstrates that the attenuation effect on the 
observed kinematics becomes stronger for larger values of 
ojoc, which agrees with the prediction that more scattering 
events correspond to a larger effect on the observed kine- 
matics. Notice that for the reasonably large gradients in lo 
(more than 30 per cent), the differences between the vari- 
ous kinematical profiles are relatively small, and definitely 
within the observational errors. 



6 DISCUSSION 

6.1 Dark matter halos around elliptical galaxies 

In the outer regions of the galaxy, the observed kinematics 
are strongly affected by photons emitted in the central re- 
gions of the galaxy, that leave the galaxy after one or more 
scattering events along lines of sight with a large projected 
radius. Because these photons are emitted by stars that gen- 
erally have larger velocities than the typical line-of-sight ve- 
locities appropriate along these lines of sight, they bias the 
LOSVDs towards larger velocities, and cause high-velocity 
wings. As a result, the projected velocity dispersion and the 
shape parameters ki and /14 increase significantly at large 
projected radii. 

These results are particularly important for the inter- 
pretation of the stellar kinematical evidence for dark matter 
halos around elliptical galaxies. For disc galaxies, the obser- 
vational evidence for the existence of dark matter halos is 
convincing: the Hi rotation curves that remain flat or even 
rising out to very large radii, provide a clear proof of their 
existence (Freeman 1970; Faber & Gallagher 1979). For el- 
liptical galaxies, this important tracer can generally not be 
used. A number of early-type galaxies, most of them classi- 
fied as SOs, have neutral or ionized gas discs which can be 
used to estimate their dark matter content (Bertola et al. 
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1993; Franx et al. 1994), but these galaxies are exceptional 
cases and may not be representative for the general class 
of elliptical galaxies. The most convincing evidence for the 
existence of dark matter halos comes from measurements of 
the density and temperature of their hot X-ray emitting at- 
mospheres (Forman, Jones & Tucker 1985; Matsushita et al. 
1998; Loewenstein & White 1999), and from gravitational 
lensing (Griffiths et al. 1996; Keeton, Kochanek & Falco 
1998). This evidence indicates that elliptical galaxies must 
have very large mass-to-light ratios at large radii, but, un- 
fortunately, they do not contain much information about the 
detailed structure of a dark matter halo and its coupling to 
the luminous matter. 

The most important way to trace dark matter halos 
around elliptical galaxies at small scales is by studying the 
stellar kinematics. With the present 8m class telescopes, stel- 
lar kinematics can be reliably traced out to several effective 
radii. If a dark matter halo is present, one expects the ve- 
locity dispersion profile to drop only slowly or to remain 
constant with projected radius. Such a behaviour was inter- 
preted as a signature for the presence of a dark matter halo 
(Saglia, Bertin & Stiavelli 1992; Saglia et al. 1993). However, 
a slowly decreasing or nearly constant velocity dispersion 
profile can also be due to a strong tangential anisotropy at 
large radii. This so-called mass-anisotropy degeneracy (Ger- 
hard 1993) can be broken by studying the LOSVD shape 
parameters: galaxies with a tangentially anisotropic orbital 
structure have a negative hi. The combination of a slowly 
decreasing velocity dispersion profile and a positive /u at 
large projected radii, is generally interpreted as a indica- 
tion of the presence of a dark matter halo. For a number 
of elliptical galaxies, the existence of dark halos has been 
advocated by such evidence (Rix et al. 1998; Gerhard et al. 
1998; Kronawitter et al. 2000; Gerhard et al. 2001). 

We show here that scattering off dust grains has the 
same effect on the dispersion profile as a dark halo: the dis- 
persion will decrease more slowly than expected. Moreover, 
also the /14 profiles are considerably biased towards larger 
values, such that the signature of an intrinsically tangential 
anisotropy can be weakened. Dust attenuation hence makes 
a foolproof detection of dark matter halos from stellar kine- 
matical evidence much more complicated. It is obviously 
important to investigate to which degree dust can reduce 
or eliminate the need for a dark matter halo to explain the 
observed kinematics. Our first results indicate that a dust 
component that is shallower than the stars with tv = 1, has 
the same kinematic signature as a dark matter halo that con- 
tains half of the total mass of the galaxy (Baes & Dejonghe 
2001c). An in-depth investigation of this topic, however, is 
beyond the scope of this paper, and a forthcoming paper 
will be devoted to this problem. 

6.2 The determination of the dust distribution in 
elliptical galaxies 

In the previous subsection, dust attenuation is considered 
as something troublesome or inconvenient - as an obstacle 
that prevents the observation of the true projected kinemat- 
ics, and therefore complicates the interpretation of the true 
dynamical structure of elliptical galaxies. The fact that the 
observed kinematics are seriously affected by dust attenua- 
tion, can also serve in a more positive way: it can be used to 



trace the spatial distribution of the diffuse dust component 
in elliptical galaxies. 

This is illustrated in Fig. |^, where we illustrate how 
the observed kinematics change with colour. The a p and hi 
profiles in blue bands are more strongly affected by dust at- 
tenuation than those in red or near-infrared bands, because 
the extinction efficiency decreases with wavelength (Table 
|l|). Combining this with the sensitive dependence of the ob- 
served kinematics on the dust geometry, we see that multi- 
colour kinematical profiles can help to constrain the spatial 
distribution of the dust. Indeed, if the dust is shallower than 
the stars (a small), the blue kinematics will be significantly 
affected by dust attenuation, whereas this effect will gradu- 
ally become weaker if we look towards longer wavelengths. 
If, on the other hand, the dust has a similar distribution 
as the stars, the observed kinematics are hardly affected, 
even in blue colours, such that we will see no change of the 
kinematics with wavelength. 

This method is also subject to a degeneration: differ- 
ences between the stellar kinematics at different wavelengths 
can also be due to variations in the stellar populations, i.e. 
different absorption lines can trace different stellar types 
which do not necessarily need to be in the same dynami- 
cal state. The major problem, however, is the observational 
challenge: the projected kinematics need to be measured at 
different wavelengths with a sufficient accuracy out to sev- 
eral effective radii. Whereas this is nowadays possible in the 
optical, it is probably beyond the limit of the possibilities of 
the current generation of telescopes and instruments to do 
this at near-infrared wavelengths. 

As a remark, we want to note that a similar principle 
has already been applied in another context: a number of au- 
thors have tried to constrain the dust content of disk galax- 
ies by multi-wavelength analysis of the apparent rotation 
curves. Bosma et al. (1992) demonstrated how dust absorp- 
tion can affect the observed rotation curve of edge-on spiral 
galaxies, and argued that the comparison of the optical and 
Hi rotation curves can be used as a new opacity test. Ap- 
plying this test to the edge-on spiral galaxy NGC 891, they 
found that at least the outer regions of this galaxy must be 
optically thin. In a similar way, Prada et al. (1994) mea- 
sured the ionized gas rotation curve of the inclined galaxy 
NGC 2146 at different wavelengths. They found a few dis- 
crepancies between the optical Ha and the near-infrared 
[Sill] rotation curves in the center of the galaxy, which they 
could attribute to a dust lane, and a nice agreement be- 
tween the rotation curves at large projected radii. Whence 
they concluded that NGC 2146 has to be largely transparent 
in its outer regions. 

6.3 The central velocity dispersion in dusty 
galaxies 

In Paper I, we found that the effects of absorption on the 
observed kinematics are only noticeable in the most central 
regions. For optical depths of order unity, absorption causes 
a decrease of the central dispersion with a few per cent. 
When scattering is taken into account, the effect of dust 
attenuation on the central dispersion is more complicated 
to predict. In particular, the nature (increase or decrease) 
and strength depend on the orbital structure of the galaxy. 
Although the magnitude of these effects is still modest (in 
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Figure 9. The observed kinematics in different colours as a tool to constrain the dust distribution in ellipticals. This figure shows how 
the projected velocity dispersion profile (upper panels) and the /14 profile (lower panels) as they are observed in different wavebands. 
In each panel, the wavebands shown are the U, B, V, I, J and K bands. All models correspond to an isotropic Plummer model with a 
V-band optical depth ry = 1, but the three columns correspond to three different dust geometries: a = f.25 (left), a = 3 (template 
model, middle) and a = 5 (right). 



particular compared to the effects of dust attenuation at 
large projected radii) , this has some important implications. 

A first area where the effects of dust attenuation should 
be taken into account are mass estimates. Indeed, the ob- 
served central velocity dispersion of galaxies is a parameter 
that is often adopted to quantify the total mass of a galaxy. 
For elliptical galaxies in the nearby universe, more accurate 
dynamical mass estimates are available through modelling of 
the entire observed kinematics. For other systems, however, 
in particular in the high-z universe, the central velocity dis- 
persion is often the only available kinematic property that 
can be measured with some reliability. The effects of dust at- 
tenuation can hence lead to a bias in simple mass estimates. 
Moreover, it is currently unclear how much dust is present 
in high-z galaxies - a significant amount of dust production 
can be expected at the early epochs of star formation. 

The central velocity dispersion also appears as an im- 
portant parameter in the discussion about the presence and 
masses of black holes. It is believed that a significant frac- 
tion of the galaxies, if not all, harbour a massive black hole 
in their inner regions. The most reliable determination of 
the black hole masses in quiescent galaxies is by means of 



spatially resolved kinematics. This, however, represents a 
challenge for both observers and modellers, because it re- 
quires a very high spatial resolution and detailed dynami- 
cal modelling techniques. Mainly with high-resolution HST 
data, black hole masses have been determined for a num- 
ber of nearby galaxies. It was found that they are tightly 
correlated with the central velocity dispersionn of the host 
galaxies by a relationship M, oc a a , with a ~ 4 (Gebhardt 
et al. 2000a; Ferrarese & Merritt 2000). Black hole masses 
of nearby AGNs, which can be determined by reverbera- 
tion mapping (Blandford & McKee 1982; Kaspi et al. 2000), 
seem at first sight to satisfy this relation as well, albeit with 
some more scatter (Gebhardt et al. 2000b; Ferrarese et al. 
2001). It is of course very tempting to adopt this relation 
for the determination of black hole masses in other galaxies, 
in particular at high z, where spatially resolved kinematics 

6 Gebhardt et al. (2000a) use an effective dispersion in their re- 
lation, defined as the aperture dispersion within lR e , whereas 
Ferrarese & Merritt (2000) adopt an aperture dispersion within 
an effective aperture of radius R e /S. 
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are beyond the present observational capabilities. However, 
because the lack of knowledge about the dust content of 
these galaxies, some caution is advised when applying this 
relation without taking attenuation effects into account. 

6.4 Disk heating processes in spiral galaxies 

It is generally known that the velocity dispersion of stars in 
the Galaxy increase with age (Wielen 1977). This is believed 
to be a result of the gradual heating of a initially cold disk 
due to irregularities in the gravitational potential. A num- 
ber of possible scattering agents have been proposed for this 
heating, but the two most important mechanisms that con- 
tribute to this heating are thought to be scattering off giant 
molecular clouds (Spitzer & Schwarzschild 1951, 1953) and 
spiral density waves (Carlberg & Sellwood 1985). The study 
of the velocity ellipsoid in disk galaxies provides an interest- 
ing way to distinguish between these mechanisms, because 
each of them leaves a different kinematical signature. Spiral 
density waves are inefficient in scattering the stars in the 
vertical direction, such that it will result in a small <j z /or 
ratio. On the other hand, molecular clouds tend to scatter 
stars rather isotropically, yielding an intermediate value for 
o~z/(JR. The axis ratio of the velocity ellipsoid is hence a 
very useful tool to identify the principal heating agent in 
spiral galaxies (Jenkins & Binney 1990; Merrifield, Gerssen 
& Kuijken 2001). 

Unfortunately, the determination of the shape of the 
velocity ellipsoid in spirals is not straightforward. Gerssen, 
Kuijken & Merrifield (1997, 2000) showed that, in theory, 
it is possible to constrain the shape of the velocity ellipsoid 
of intermediately inclined spiral galaxies from the observed 
dispersion profiles on the major and minor axes. 

From the results in this paper, however, we anticipate 
that dust attenuation has a strong effect on the projected 
dispersion profiles in disk galaxies, because of two reasons. 
First, disk galaxies contain large amounts of dust, typically 
several orders of magnitude larger than the average ellip- 
tical galaxy. The effects of attenuation are therefore likely 
too important to be considered as a second-order effect. Sec- 
ond, the differences in intrinsic velocity dispersion between 
bulge and disk are very large. If a photon, emitted by a 
high-velocity bulge star, propagates into the disk, and is 
scattered there such that it leaves the galaxy at a large pro- 
jected radius, it will contribute a very large line-of-sight ve- 
locity to that LOSVD. A relatively small number of bulge 
star photons can therefore already seriously contaminate the 
LOSVDs and increase the observed projected dispersion at 
these lines of sight. Moreover, both the dust content and the 
bulge-disk ratio vary along the Hubble sequence, which def- 
initely complicates a simple picture. These ideas definitely 
need a detailed investigation before major conclusions can 
be drawn on the heating processes in disk galaxies from the 
observed velocity ellipsoids. 



7 CONCLUSIONS 

The aim of this series of papers is to investigate the effects 
of a diffuse dust component on the observed kinematics of 
elliptical galaxies. We started in Paper I by investigating the 
effects of absorption only, neglecting the scattering effects. 



We found that, for realistic optical depths, these effects are 
modest, i.e. of the order of a few per cent in the central 
regions and completely negligible at larger projected radii. 

In this paper, we extended our models to include the ef- 
fects of scattering, which are usually considered as a second- 
order effect compared to absorption. The underlying thought 
is that the main effect of scattering is a reduction of the ef- 
fects of absorption. It is therefore believed that the effects of 
scattering can be modelled by considering absorption with 
a reduced effective optical depth. From this point of view, 
it is often considered not worth the effort to properly in- 
clude scattering in a proper way, thereby bypassing a costly 
radiative transfer treatment. 

If we would have adopted this prescription, there would 
have been no reason to include scattering in our models. In- 
deed, because the effects of absorption only on the observed 
kinematics are already fairly small, the extra second-order 
effect of scattering should be completely negligible. How- 
ever, a number of authors, foremost Witt et al. (1992), have 
shown that the effects of scattering are important, even for 
small optical depths, and that any way of neglecting or ap- 
proximating them can lead to serious errors. We confirmed 
this in a previous study on the effects of dust attenuation 
on disk galaxies (Baes & Dejonghe 2001b). 

In this paper, we have clearly demonstrated that, con- 
cerning the observed kinematics in elliptical galaxies, scat- 
tering can not be considered as a second-order effect to ab- 
sorption. On the contrary, we find that the effects of dust 
attenuation are much more complicated and fascinating if 
scattering is included in the modelling. The way the kine- 
matics of elliptical galaxies are affected can differ drastically, 
depending on which line of sight is considered, on the star- 
dust geometry and on the internal orbital structure of the 
galaxy. The most striking effect is the serious increase of 
both the velocity dispersion and the LOSVD shape param- 
eters at large projected radii, due to photons from high- 
velocity stars scattered into the line of sight. This effect, 
which complicates the interpretation of the stellar kinemat- 
ical evidence for dark matter halos around elliptical galaxies, 
has absolutely no counterpart when only absorption is taken 
into account. Results as these, which may seem unexpected 
at first sight, can be easily understood, once the idea of scat- 
tering as a second-order effect of absorption is abandoned. 
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APPENDIX A: THE SPATIAL LOSVD FOR 
THE PLUMMER MODEL 

Al The spatial LOSVD 

The dynamical structure of a certain stellar population of 
a galaxy is completely determined by its phase space dis- 
tribution function F(r,v). Whereas this function gives us, 
at a certain position in the galaxy, the three-dimensional 
distribution of the velocities, it is interesting from an obser- 
vational point of view, interesting to know the (marginal) 
distribution of velocities into an arbitrary direction k, for 
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example the direction towards the observer. To calculate 
this distribution, we construct a new cartesian reference sys- 
tem (ei^eij.eii), such that fc = en. The spatial LOSVD 
cj>(r,k, v») is defined through the relation 



(Al) 



Z(r)<t>{r, k,v\\) =C II F(r,v) dv± ± dv±_ 



The emissivity is taken into this definition in order to have 
the normalization 



fr,fe,?;ii)d«ii 



C 

£(r) 



F(r, v) dv 



(A2) 



The spatial LOSVD hence represents the probability for a 
star at a position r, to have a velocity component vu in the 
direction fc. The quantity £(r) 4>{r, fc, v« ) then represents the 
stellar emissivity at r of stars which have a velocity vv, in 
the direction fc. 

For the calculation of the spatial LOSVD, it is impor- 
tant to know that a distribution function will not depend 
on the six phase space coordinates (r, v) independently, but 
only through a number of integrals of motion. In a spherical 
galaxy, we can write the distribution function generally as 
F(E,L), where E and L are the binding energy and angu- 
lar momentum integrals respectively (Binney & Tremaine 
1987). Because the distribution function only depends on Vg 
and through the combination v% = vj+v 2 ,, we can choose 
our reference system such that e± 2 — eg. If we denote the 
angle between fc and e r as n, we then find 



E = 1>(r) - \vl, 



1 2 

9 V ±2 



1 2 
2 V \\ 



2/2 2 2 2-2 

r (v± 1 cos 77 + vj_ 2 + sin rj ■ 



(A3) 

«H sin 277) . (A4) 



Substitution of these expressions into the distribution func- 
tion allows i n pr inciple to calculate the spatial LOSVD using 
expression (Al). 



A2 The Plummer model 

A major advantage of the Plummer model is it allows the 
construction of a completely analytical one-parameter fam- 
ily of self-consistent dynamical models (Dejonghe 1987). The 
different models are characterized by a parameter q, re- 
stricted by q ^ 2. This parameter is determines the or- 
bital structure of the galaxy: negative values of q corre- 
spond to a tangential anisotropy, positive values to a radial 
anisotropy, whereas for q = the model is isotropic. By 
means of Laplace-Mellin transforms the corresponding dis- 
tribution function F(E, L) can be calculated analytically for 
all values of q (Dejonghe 1986). For general values of q, the 
distribution function can be expressed in terms of hyper- 
geometric functions. For even integer values of q however, 
the distribution function becomes much simpler, which will 
allow us to calculate the spatial LOSVD explicitly. 



A2.1 The isotropic model q — 

For the isotropic dynamical models, the distribution funci- 
ton is only function of the binding energy. For the Plummer 
model, it is a simple power law of E, 



F(E) 



_3 1_ 

7n 3 c 3 it 3 



•IE 



7/2 



(A5) 



where u = GM/c represents a cha ract eristic velocity. Sub- 
stituting the distribution function dA5f) in to the definition 
(|Al|), we obtain, using the expression ( pV3[ ) , 



(r,k,vn) = 



63tt y^r) 



1 



9/2 



(A6) 



Because the distribution function of isotropic galaxies is 
completely symmetric in the three velocity components, 
the spatial LOSVD is independent of the direction fc. It is 
straightforward to check that the normalization condition 
(A2) is satisfied. 



A 2. 2 The radial model q = 2 

For q — 2, the most radial model in the Plummer family, 
the distribution function vanishes for IE ^ L 2 /c 2 , whereas 
reduces to a simple form for 2E L 2 /c 2 , 



T 2 



3/2 



(A7) 



To calculate the spatia l LO SVD, one has to substitute the 
expressions ( |A3| ) and (A4) into this distribution function, 
and integrate it with respect to v± 1 and v± 2 , where the 
integration surface in determined by the condition 2E ^ 
L 2 /c 2 . After some algebra, one finds that the LOSVD is 
given by 

4>{r, fc, vn 



16 



5tt ^/2Aip{r) 



2Aip{r) 



5/ 2 



(A8) 



for |?;|| I ^ \/2A r ip(r), with 

22 2 
r cos rj + c 
A. ~ — 1. 



(A9) 



The normalization condition (A2) can be easily checked 



A2.3 Tangential models with q 



-2m 



For general negative values of q, the Plummer distribution 
function can be written in terms of hypergeometric series. 
If q is a negative even integer, q = —2m, this series breaks 
down after a finite number of terms, and the distribution 
function then becomes a finite series, where each term is a 
power law of E and L (Dejonghe 1986). The spatial LOSVD 
can be calculated for each of these distribution functions 
following the same recipe as for the radial q — 2 model. An 
alternative, elegant way to calculate the spatial LOSVDs for 
these models using Laplace transforms is presented by De 
Rijcke (2000). After lengthy calculations, one finds eventu- 
ally 



{r,k,vu 



1 r(6 + 2m) 1 
0FT(4i+2m) ^2~W) 



2tp(r) 



3=0 



, 2-2 \ 2j 

j 1 r sin r\ x 



(JO 2 



1 - 



2ip{r) 



x 2 J\ -j,5 + 2m;i; 



2ip(r 



(A10) 
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Figure Bl. The probability density </>(C) f° r C = \\hi ~~ 
corresponding to the Henyey-Greenstein phase function. Curves 
are shown for different values of the asymmetry parameter g, 
ranging between g = and g = 0.8, in steps of 0.1. The black 
lines indicates the probability density for g = 0.6, appropriate for 
the optical regime. 



where (a)j represents the Pochammer symbol 
(o)j = o(o + l)(a + 2).--(o + j-l). 



(All) 



Setting q = m = in the expression ( A1C ), it reduces to the 
expression (A6) we found for the isotropic Plummer galaxy. 
To check the normalization of the spatial LOSVDs we have 
to evaluate the integral 



l(m,j) 



f 

Jo 



(1 



2)1+2— i 2Fl (-j, 2m + 5;±;x 2 )dx, 



(A12) 



for all j ^ m. By means of formula (7.512.4) of Gradshteyn 
& Ryzhik (1965) we find that this integral vanishes for j > 0, 
whereas for j = we obtain 

r(^+2m) 
2 T(6 + 2m) ' 



l(m,Q) 



(A13) 



Only the first term in the summation in equation ( A10) will 
hence contribute to the integral over v\\ , and it is obvious 
that the final result will equal 1. 



APPENDIX B: ANISOTROPIC SCATTERING 
AND THE IMPORTANCE OF THE DUST 



GRAIN VELOCITY TERMS IN EQUATION (|10|) 

Each term in the summation in expression ([To]) is a inner 
product of the velocity vector Vd f and the relative direction 
vector hi — fcj_i. It is obvious that the contribution of one 
such term does not depend only on the dust grain veloc- 
ity u t j i , but also on the scattering angle. The distribution of 
scattering angles will therefore also play a role in the im- 
portance of the dust grain velocity terms in equation (|lo|). 
If, for example, the propagation direction of a photon is on 
average only slightly changed during a scattering event, we 
find that hi ~ fe»_i, which reduces the magnitude of the 
dust grains velocity terms. 

We can illustrate this nicely by calculating the statis- 
tical distribution of the modulus of the vector hi — ki-i, 



Table Bl. The mean, the median and the maximum value for the 
probability density ip(C) corresponding to the Henyey-Greenstein 
phase function, for different values of the anisotropy parameter 
9- 
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Cmcd 


Cpcak 


0.0 


1.33 


1.41 


2.00 


0.1 


1.25 


1.30 


2.00 


0.2 


1.16 


1.19 


1.26 


0.3 


1.07 


1.06 


0.90 


0.1 


0.98 


0.93 


0.67 


0.5 


0.87 


0.79 


0.50 


0.6 


0.75 


0.64 


0.37 


0.7 


0.63 


0.49 


0.25 


0.8 


0.48 


0.33 


0.16 


0.9 


0.29 


0.17 


0.07 


1.0 


0.00 


0.00 


0.00 



i.e. we will calculate the probability density <p(() for C = 
\\ki — fei_i||. Since the directions fc;-i and hi represent unit 
vectors, elementary trigonometry learns that 



c = - 



(Bl) 



where a is the scattering angle. The scattering angle is dis- 
tributed according to the phase function $(cos a), such that 



1 - 



C 



(B2) 



with C taking values between and 2. If the scattering is 
isotropic, the phase function has the simple form "I>(cos a) = 
1, such that <^(C) = C/2- Larger values for ( will hence 
be pref erred. If we adopt Henyey-Greenstein scattering (see 
Section 4.2.2), we can calculate the corresponding values for 



different degrees of anisotropy. In Fig. Bl, the probability 
density <fi(Q corresponding to the Henyey-Greenstein phase 
function ( [l7[ ) is shown for different values of anisotropy pa- 
rameters g. In the limit g = 0, this function of course reduces 
to the isotropic equivalent (p(() = C/2- If the anisotropy pa- 
rameter increases, however, smaller values for £ will gradu- 
ally become dominant. We can quantify this trend by calcu- 
lating the mean value Cmean, the median ( mc d or the value 
Cpcak for which ip(C) reaches its maximum value. For the 
isotropic phase function we find C m0 an = f , Cmcd = %/2 and 
Cpcak = 2, whereas for the Henyey-Greenstein phase func- 
tion, we obtain 



Cn 



1-3 



l±£ ]n (l+^ 



- 1 



= VV - 3.g + 2 
1-.9 



Cmcd 

Cpeak = min | 2 



(B3) 
(B4) 
(B5) 



In Table Bl we tabulate Cmean, Cmcd and Cpeak as a func- 
tion of the asymmetry parameter. This table confirms that 
smaller values of C quickly start to dominate as g grows 
larger. In particular, for the anisotropy parameters appropri- 
ate in the optical regime, g ~ 0.6 (see Table ^), the density 
<p(C) clearly favours fairly small values for C (black curve in 
Fig. Bl). As a consequence, ||fei — fct— 1|| will be fairly small in 
the majority of the scattering events. The anisotropic nature 
of scattering hence contributes to reducing the importance 
of the dust grain velocity terms in equation (JlCl) . 



This paper has been typeset from a TgX/ WTftX. file prepared 
by the author. 
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